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1. INTRODUCTION 

Simultaneous wireless information and power transfer (SWIPT) networks are classified into two 
categories: direct and cooperative relaying system. The relays are employed with several protocols which 
help to guard them against extreme signal fading, i.e., amplify-and-forward (AF), and decode-and-forward 
(DF). Moreover, these relaying networks operate in full-duplex (FD), half-duplex (HD) or even hybrid 
FD/HD mode, allowing the system to switch opportunistically between the two modes, which outperforms 
both the standalone modes, [1]-[3]. To further improve the network performance, theories about opinion 
dynamics [4] and block chain for security issues [5] could be promising drivers. 

Besides, network designers have to consider as well the fading modelling techniques for the 
wirelesschannel to find the best fit for the simultaneous wireless information and power transfer (SWIPT) 
systems in practice [6]. Log-normal fading channel, indeed, showsit strength in characterizing indoor 
shadowing effect of human body movements, as well as obstacle mobilitiesand walls of buildings [7]-[12]. 
Thus, it is an appropriate model for modelling shadow effect for internet of things (IoTs) and smart city 
applications despite its drawback in modelling the small-scale fading in low-power sensor network. In the 
context of the aforementioned indoor scenarios, the utilization of cooperativerelaying system becomes 
inevitably a necessity. It is beneficial not only for increasing the channel capacity but also for spreading the 
network coverage range while opposing the shadowing effect [13]. Nevertheless,the higher number of relays 
means the higher level of inter-relay interference problem, which requires relay selection (RS) schemes to 
solve [14]-[16]. Moreover, it is obvious that battery plays an essential role in improving the system 
performance of cooperative battery-operated relaying networks. In such network scenarios, in [17]-[20] 
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investigated the dynamic behaviours of finite-capacity battery model and its topology utilizing the property 
of the Markov model. In particular, the batteries at the relays have some initial energy and will harvest 
additional energy during the transmission process. Additionally, the battery status is considered when 
choosing the best relay. 

All the aforementioned technologies, although being applied solely for wireless networks, can be 
combined with household low-power energy harvesting (LoPEH) wireless sensor networks (WSN) to bring 
about their advantages for smart homes and smart cities development. A limited number of studies have 
proved that utilization of relays can effectively eliminate the LoPEH WSN’s problems of high noise and 
signal attenuation [21]-[23]. This opens up a bright future for 5G development especially for indoor 
applications, multi-building LoPEH WSN, whereas there are underground constructions, metal and concrete 
walls [24]-[25]. 

With inspiration taken from the above works, this paper conducts the performance analysis of relaying 
SWIPT systems in terms of the overall outage event (OE) with a relay selection (RS) scheme so-called optimal 
relay selection (OPRS) for multi-relay cooperative wireless networks in FD-AF protocol. The LoPEH WSN 
uses battery-operated relays and the shadowing effects were simulated with log-normal fading. Thanks to the 
Markov property of the energy buffer status, the stochastic property of the harvested energy flow can be 
characterized with a proposed active-passive model. The rest of this paper is arranged as follows: i) section 2 
presents the system model, ii) section 3 analytically derives the performance indicator of the proposed the 
Markov battery model, and iii) section 4 subsequently presents the numerical results. Finally, this paper is 
summarized in conclusions. 


2. SYSTEM MODEL 

A cluster network is illustrated in Figure 1 which includes a source (S), a destination (D), 
and a cluster (C) of in-between K cooperative relays (R;) (1 <i < N = K) installed in a LOPEH WSN. 
The assumed system is under coverage extension scenarios, whereas the (S) and (D) nodes are considerably 
distant from each other, which makes the direct link between them so attenuated that the communication 
cannot be realized without relays [26]. Besides, cooperative relaying networks can be utilized as well for 
situations where the surrounding environment causes deep shadowing effects on the direct link with physical 
obstacles’ mobility [27]. Furthermore, it is assumed that the carrier and symbol synchronization are ideal, and 
each terminal knows the channel state information (CSI) in advance. Information flows from (S) node to (D) 
node given that (S) is powered by a conventional stable source being P, and (R) node is powered by a battery 
Pp implemented with the energy harvesting (EH) module. In addition, there is an additive white Gaussian noise 
(AWGN), nj, j E {r, d} with zero mean and variance No respectively at i — th (R) and (D) nodes. 


R À Makov Prisma] 
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ma S-D channel 


Figure 1. System models 


Besides, Figure 1 schematically presents the distances between (S) — (R), (R) — (D), and (S) — (D) 
nodes which are denoted by dsr, dya, and d,q with corresponding channel coefficients lsr, l-a, and loa- 
Similar to the configurations for studying the log-normal fading channels in [28], the (S) and (D) nodes are 
each equipped with an antenna. Every (R) node is implemented with AF protocol and equipped with two 
antennas (one for receiving and the other for transmitting purpose). As aforementioned, this FD operation 
inevitably introduces the loop interference channel to the system indicated by l,„y. Particularly, the 
communication time is split into slots given that within each time slot, the i — th relay R; , (R; E N) is 
selected according to some relay selection (RS) scheme to aid the (S) node’s transmission. In a signal block, 
s(t) denotes the narrow-band transmit signal at (S), (R) has zero mean, and E[|s(t)|?] = 1 with E denoting 
the statistical mean operation. 
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ee there are assumed independent and identically distributed log-normal random variables 


(RVs) 12, ed and a , respectively assigned with parameters LN (201 40? ), LN (20, 44%. 4); and 


Isr 
LN bap a ADI, x The loop interference channel l2, under log-normal distribution is assigned with 


parameter LN (2,,.,. 407). This parameter plays an essential role in characterizing the loop interference 


strength, and thus, determines the overall performance of the FD relaying network. 

Remarkably, (S) — (R) and (R) — (D) links are under the influence of quasi-static block fading. 
This means that the channels are constant at time block T, then vary from block to block as they are 
independently and identically distributed (i.i.d.) in log-normal fashion. The probability density function 
(PDF) and the cumulative distribution function (CDF) of the RV X in log-normal distribution are respectively 


calculated by: 
2 
(mas In(z)-20x) 
0/In(40) , ee 


2 
Z /8TOAXY 





I n(z)-20 x 


Fy(z) =1-Q (ae) and fy(z) = 





t2 
is a scaling constant and Q(-) is the Gaussian Q-function, Q(x) = ie es Aar 


Inwhich § = T 


T 
2.1. Direct link without cooperative protocol 

As aforementioned, the outage event (OE) is an indicator of the system performance. This is the 
probability that the instantaneous capacity falls below a given threshold bits per channel use (BPCU) being 
Ro, and Pr|Cs a < Ro]. Hence, by utilizing the CDF of the log-normally distributed RV |l a ie the OE of the 
direct link can be expressed as: 


OF.4 = Pr (isal < (2¥-1) zei) =1-@ (ja) a) 


2054 





Where a = rah and yọ = 2%0/W _ 1, m is the path loss exponent, Ry is the data transmission rate, W is 
s' 


the frequency bandwidth. 


2.2. Relay-aided cooperative protocol 

Definition 1. The relay-aided cooperative protocol is utilized to further improve the ergodic 
capacity, hence the LoPEH WSN performance. If the direct link is considerably attenuated with deep fading, 
an alternative relay link will be selected among Ks to conduct the transmission instead, on condition that it 
possesses enough energy for the task. 

Accordingly, the time frame T for completing one transmission cycle from (S) to (D) is split into 
three slots according to the time-switching relaying (TSR) protocol as depicted in [16]. The first time slot 
being tT (t is EH time factor, 0 < tT < 1) is dedicated to EH purpose whereas (R;) can harvest the energy 
from the signal broadcasted from the (S) node. The rest is halved with (1 — t)T /2 utilized for (S) — (R;) and 
(R;) — (D) information transmission. 
neT Ps|ls.r,| 


m $ 
ds ri 


Within the first time slot, based on the energy harvested at the i,, (R) node, Ey = 


the relay’s transmit power can be formulated as: 


2 
p= Ey = nrPs|lsr,| 
Ri a-vrT (1-1) a, 





(2) 


Where 0 <n < 1 is the FH efficiency specified mainly by the characteristics of the circuitry. 
For the second time slot being the information transmission phase, the received signal at the itn (R) 
can be calculated by: 


y(t) = ve le lsr;S(t) + L,rr(t) + Ny (3) 


Where the information signal s(t) is normalized as E[|s(t)|?] = 1, and the loop interference r(t) of the FD 
relays satisfies E[|r(t)|*] = Pr- 
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It is noteworthy that a relay in a FD multi-relay cooperative system can distinguish its signal with 
others’ [29]. Therefore, it can apply interference cancellation for its own loop interference reduction. 
The post-cancellation signal at the relay, (3) can be expressed as: 


yr (t) = faa ler S(t) + Ur P(t) + ry (4) 


Whereas ae is the residual loop interference remaining after imperfect interference cancellation process, and 
E{lA()I?} = Pr, 

In FD-AF system, after base-band processing at (R;) following (4), the signal is amplified and 
forwarded to (D). Accordingly, the (D) node receives the signal given by (5): 


PsPR. PR; PR; 
Ya (t) = =e leritria Gs(t) ae Fae Lyd Lr Gr(t) T m lra Gn, + Nq (5) 
dsr, 4rd drid dria 


Where the relay gain G normalizes the instantaneous received power while permitting the relay transmission 
with maximum power as: 


1 1 
P, > œ, G = z% 


Ps 2 — 2 Ps 2 r2 
[ath +|isri| PR;+No [abled +isri| PR; 








For such case, after substituting G into (4) some basic algebraic manipulations are performed to 
obtain the end-to-end signal-to-noise ratio (SNR) of the i,;, relay at (D), as in [2]. 


PSYS,R; p 
_  PRiYRRį bie a 6 
Ys,Rid = PSYS,R; Tp a ( ) 
PRiYRRįi RiYR;,D 
Where: 
2 2 2 
= -m = =m = 
YS,Ri AI [Hox | dsr» YR;,D AEN [l-al drd and YR,Ri TE Lea 


Accordingly, it should be noted that the instantaneous capacity of the FD-AF-TSR system is obtainable from: 


Cyr,a = (1-1) log,(1 + Ysa) (7) 


Here, RS schemes for the relay-aided cooperative LoPEH wireless sensor networks, namely optimal 
relay selection (OPRS) is established with regard to (6). Thus, the best relay k,, is activated when 
k = arg maxfYs r,a} The best selected relay is responsible for establishing the communication and subjects 

U 


to the performance analysis. 


3. PERFORMANCE ANALYSIS 
3.1. Energy storage modeling at relay 

Thanks to the stationary stochastic process stated in [30], ones can denote the minimum energy 
required for Rẹ activation and the energy which is harvested from the environment during the period of the 
kın signal block as en,. Particularly, the very first focus is paid on Rẹ without energy storage working solely 
on enx. The relay runs out of energy if eng < PR. The probability of such events is so-called energy-exhausted 


probability and expressed as follows: 
P 
OE, = Prien, < Pp, | = J, ™ fen, (x)ax (8) 


Where fen, (x) signifies the PDF of the stationary stochastic process being eng. 

On the other hand, for most of the applications in practice, the Rx is installed with a battery whose 
capacity is assumed to be enmax. The energy consumption for the ktp signal block transmission is denoted as 
OE,. With regard to the relay-aided cooperative protocol scenario described in Definition 1 , ones can 
express the OE, with the stationary random function as follows. 
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OE, a/K if Pr(oup = Pry) 
Pr(our)=5\ 1-—0E,a/K, if Pr(ou, = 0) (9) 
0 otherwise 


Where factor K indicates that every (R) node is activated with an equal probability. 

It is noted that en; and ou, are two stationary and independent variables. To describe the energy 
buffer status following the Markov stochastic process, S% is utilized, and it indicates as well the initial energy 
amount stored at (R,) when the k,, signal block transmission starts. Accordingly, the stationary PDF of Sk 
being fs„(x) is obtained as a combination of (8), (9) and the original PDF of eng. Theory-wise, f,, (x) is 
sufficient for describing all the stochastic characteristics of the EH module. Nonetheless, this discretized 
model has a large number of possible energy buffer statuses. Therefore, a simplified model is proposed to 
limit the number of statuses to two processing steps in the Lemma 1. 

— Lemma 1. To decide whether the direct link or the alternative relaying link is going to be utilized for 
data transmission, S is compared with the possible energy consumption at (Rẹ) as: i) If sy = Pry» 
the kıp relay is activated and relay the signal for (S)-(D) link on condition that the direct (S) — (D) link 
is too weak to perform its task and ii) If sz < Pry» the ky, relay remains inactivated and consumes no 
energy. Communication is realized via direct link. 

Practically, these two statuses can be utilized for characterising the flow of harvested energy. This is 
named 1-0 (on/off) model with denotation s;, please refer to Figure 1. Thanks to the stationary PDF of Sp 
denoted as f;, (x), the stationary PDF of są can be derived as follows: 


(Seg fa Gdx, if Pr(sj) =1 
Pr(s))=4 p“ (10) 
i ei dx, otherwise 
To conclude, from the systematic point of view, a 1/0 (on/off status) model for the harvested energy 
flow is established and described in Theorem 1. For more details of how the original PDF of the harvested 
energy flow is derived and how the fen, (x) is utilized for calculating OE, with low computation complexity, 
please refer to [31]. It should be noted that this function is different from the original PDF being fen, (x) of 
the harvested energy in (9). 


3.2. Overall outage event 

It is impossible to successfully establish the relay-aided link over the Rx all the time due to the 
effect of either deep fading or energy exhausting owing to the fluctuating renewable resources. Such an event 
is so-called overall OE, and is important for performance analysis. To formulate the OE event, we describe 

first the related premilinaries in Theorem 1: 

— Theorem. 1 The relay-aided link in cooperative relaying LoPEH WSN is activated on condition that the 
direct link is insufficient for information transmission and the available energy stored at the relay node 
is greater than Pp,. Furthermore, there is an EH module employed in Rg, which is characterized with 
parameters Pp, and OE,. Thanks to this module, the stochastic property of the harvested energy can be 
captured with no loss. Given that at R,, the available energy is sufficient for signal transmission, 
the overall OE of the FD-AF-TSR system in the proposed relay-aided cooperative LoOPEH WSNs can be 
expressed as: 


OE, = OE sa X Tk ,[OE, + (1 — OE) x OE, r,a] (11) 


By considering (7), the definition of the OE for S — R, — D link in case of the OPRS scheme, 
denoted as OF, 2, 4, can be expressed in the below united manner. 


PsYS.Ry PEY 
Pag VRrg. RKYRKD 
= kik 
OE, rpa = Pry oat <y (12) 


+PRkYRęD+1 





k 
PRKYR,Rę 
Where yı = 2%0/G-)W _ 1, and the selected R, is different for the RS schemes. By utilizing the simple order 


statistic result obtained from (11) and (12), it is possible to express the overall OE of the proposed relay-aided 
cooperative transmission protocol under the OPRS scheme as follows. 
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OE (v1) = |1 - o (PS *#)| x [0E + A - OF) x OF aa] (13) 


Where: 


g 2 
oy (EIn@)-201,.,.) 
Bf Lo fein) -onr + oa) (ae) 
OF r,a) = 1- —— l zo| m he e ie i 
0 


| BTN, v2(,,,. + 1,4) 


N; Yı 1 NT 
pai G J as 
aenn 6 ee E 


, 








— Proof: here in this proof, the process to obtain the OE,, expression in (11). According to (8) from which 
R,, is selected, it is possible to rewrite the OE of the FD-AF-TSR LoPEH WSNs in (12), under OPRS 





scheme as: 
PsYS.Ry y 
PRLYR R, PK RED 
OE, R a= Pr Wee k < Yı (14) 
hee “k +PR,YR,Dt1 
PRLYRR k "k 


Noted that 6 = nt/(1—T), X = Ys, RkYry p» and Y = Yr ry: By substituting (2) into (14) the OE of the 
FD-AF-TSR LoPEH WSNs can be rewritten as: 


va(5t¥) 
OE sna = Prix < — — (15) 


Ps—Pshy1Y 
Given that X is positive, the probability in (15) can be rewritten as 
1 
¥i1(z+Y 
Pri(x> vG) =1,Y> adt 
Ps—Psbyi¥ 1 


1 
yıl z+Y 
pr(x < Tle") ve 

Ps—Psoy1Y oy1 


OE; Rid = (16) 





Besides, the complementary CDF of X, as a product of two log-normally distributed RVs, is defined as: 


1 
yi(a+x 
sola) 2(o1,,+01,4) 


Ps—Psoy1x | 


F(x) =Q a7) 


VZ(Ais rta) 


Using the PDF of log-normally distributed RV Y is obtained by: 


frx) = ——— exp (- kreni ) (18) 


2 
x nO, 801, 


Accordingly, the OE in (16) can be calculated as: 


1 1 
E (ia) 
OE sna =1- Jo" fy) Fx (aS) dx (19) 


Eventually, by substituting (17) and (18) into (19) in combination with (1). The (13) can be solved. This is 
end proof. 


4. RESULTS AND DISCUSSION 

This section presents the numerical results of Monte Carlo simulations of the overall OE of the 
relay-assisted cooperative wireless power-line communication network (PCLN) powered by a stable energy 
supplier, under the consider optimal relay selection scheme. The system parameters for the simulations are 
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as follows. The energy-exhausted event at Rẹ, OF, = 107t. Frequency bandwidth, W = 2(W). 
Transmission rate threshold, Rọ = 2 (bps/Hz). Traditional stabilized power source, P, = 5 (dB). Overall 
additive white Gaussian noises (AWGNSs), Nọ = 1. EH efficiency, 7 = 1. EH time fraction, t = 0.2. 
Path-loss exponent, m = 2. Distance from (S)—(R,), dsr =5 (m); (Ry) — (D), d,a = 5 (m), and 
(S) — (D), dsa =10(m) with corresponding channel mean Nir = 2a = Lisa = 4 (dB), and channel 
variance W, = Opa = Ola = 3 (dB). The loop interference (Rẹ) channel mean 2)... = 2 (dB) and its 
corresponding channel variance w),. = 3 (dB). It should be noted that all the channel means and variances 


are log-normally distributed. 

Figure 2 plots the overall OE of the cooperative relay-aided LoPEH WSNs versus the SNR. 
Three values of energy-exhausted probability OE, of best RS being 0, 0.5 and 1 are illustrated. It should be 
emphasized that the curve corresponds to the traditional WSNs with direct link at OF, = 1, and the relay-aided 
cooperative LoPEH WSNs with a constant power source at OF, = 0. At OF, = 0.5, the LoPEH can utilize 
both the direct link and the relay link (red, continuous curve) or only relay link (blue, continuous curve) for 
signal transmission. The overall OE curves decrease by an order of magnitude when OE% increases from 0 to 1. 
It is obvious that the relay-aided LoPEH WSN curve at OE, = 0 is the lowest within the consider SNR range, 
which indicates that it delivers the best system performance. In practice, a well-designed battery-operated 
system can have relatively small OE,. Therefore, from theory and simulation results, it is reasonable to say 
that the installation of EH relays is significantly beneficial for the WSNs. 

Furthermore, as SNR approaches 30 (dB), all the curves obtained from the OPRS scheme tend to 
converge to the lowest error floor providing zero diversity gain. The simulation results agree well with the 
theoretical analysis ensuring that the proposed scheme can be utilized further for comparison of system 
performance among other RS schemes to choose the most effective one. In addition, Figure 3 depicts the OE 
of the relay-aided cooperative LoPEH WSNs versus the SNR with different numbers of relay nodes, 
K = 1,3 and 5, in the OPRS scheme. Besides, it should be noted that all cases are plotted with OE, = 1071. 

As aforementioned, the lower the curve, the better the system performance. It can be observed that 
as K increases from | to 5, the system performance improves correspondingly. It worth to mention that in 
practice, relays embedded with EH module are relatively affordable. Hence, raising the number of on-field 
relay nodes is considered effective not only in decreasing the overall OE but also in cost manner comparing 
to other methods. As aforementioned, because the trends of the simulation results agree well with the 
theoretical analysis, they can be utilized for comparison studies between different RS schemes to choose the 
most effective one for system performance. 
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Figure 2. The overall OE versus SNR with Figure 3. The overall OE versus SNR with 
different energy-exhausted probability the number of available relay nodes K 


Figure 4 plots the overall OE behaviour under the influence of the EH parameter t in the OPRS 
scheme. The theoretical analysis and simulation results are calculated with loop interference channel variance 
DQ, = 2 (dB) and 4 (dB), in combination with the fixed system parameter P, = 5 (dB). The fact that the 
simulation results agree well with the analytical ones proves the correctness of the analysis. Additionally, 
it can be observed that as t increases to 0.15 and 0.2 respectively for DQ,» = 2 (dB) and 4 (dB), the OE 
sharply decreases. This results in a notable increase in system performance. Moreover, it is worth to 
emphasize that as T value is equal to 0 or approaching 1, corresponding to the cases when the harvested 
energy amount gets too small or too large, the OE becomes extreme, because there is no resource that 
remains for data transmission. Besides, it should be noted that the system performance can be improved by 
unequal channel allocation following the relative channel distribution method. 
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Figure 4. Compare the overall OE versus tT being the time switch 


5. CONCLUSION 

To summarize, the paper investigated the OE behaviour of a FD-AF RS scheme so-called OPRS, 
which provide zero diversity over log-normal fading channels in WSNs. Thanks to the characterized 
harvested energy flow, the two proposed cooperative relaying protocols can significantly reduce the OE 
leading to the considerable improvement of the LoPEH WSN performance. Additionally, the overall OE 
expression of the OPRS scheme was derived. Besides, as the overall OE plots versus the SNR and the EH 
time-switching in theory and in simulation results agree well with each other, the derived expressions can be 
utilized as references for comparison studies among other different RS schemes. Future papers will focus on 
analyzing the overall OE performance of EH relay-aided cooperative LoPEH WSNs under impact of 
different fading channel models along with their pros and cons. Besides, exploitation of multiple antennas 
and multiple relays distributed according to Poisson point process can be of interest to further optimized the 
existing LOPEH WSNs in every household for applications of smart homes and smart cities. 
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